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Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the 
co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body 
simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. With our 
suite of over 200 merger simulations, we identify systematic trends in the behavior of recoiling 
BHs. Our main results are as follows. (1) While BHs kicked at nearly the central escape speed 
(Vesc) are essentially "lost" to the galaxy, in gas rich mergers, BHs kicked with up to ~ 0.7 v esc m ay 
be confined to the central few kpc of the galaxy. (2) The inflow of cold gas during a gas-rich major 
merger may cause a rapid increase in central escape speed; in such cases recoil trajectories will 
£f~) depend on the timing of the BH merger relative to the change in v eS c- (3) Recoil events generally 

£f) reduce the lifetimes of bright active galactic nuclei (AGN) but may actually extend AGN lifetimes 

at lower luminosities. (4) Recoiling AGN may be observable via kinematic offsets (vbh > 500 
t-H km s _1 ) or spatial offsets (Rbh > 1 kpc) for lifetimes of up to ~ 10 - 100 Myr. (5) Rapidly- 

recoiling BHs may be up to ~ 5 times less massive than their stationary counterparts. These 
mass deficits lower the normalization of the Mbh — 0* relation and contribute to both intrinsic 
and overall scatter. (6) Finally, the displacement of AGN feedback by a recoil event causes higher 
central star formation rates in the merger remnant, thereby extending the starburst phase of the 
merger and creating a denser, more massive stellar cusp. 
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1. Introduction 

It has recently been demonstrated that "central" supermassive black holes (SMBHs) may in 
fact spend substantial time in motion and offset from the galactic center, owing to gravitational- 
wave (GW) recoil. Numerical relativity simulations of binary black hole (BH) mergers have re- 
vealed that GW recoil may impart a kick to a merg ed BH of up to 4000 km s" 1 [1]. Although 
the actual distribution of GW kick velocities is not known, estimates using optimistic assumptions 
indicate that more than a third of major galaxy mergers could result in kicks greater than 1000 km 
s~ ! [1-4]. This fraction would be smaller if the progenitor BH spins are not large, or if the spins 
could be partially aligned prior to merger via gas processes or GR precession [5-7]. 

Recoiling BHs may be observable as AGN that are either spatially or kinematically offset from 
their host galaxies [8- 13, e.g.]. To date, several offset AGN candidates have been discovered, 
though none have yet been confirmed [14-18]. We also examine possible indirect consequences 
of recoil. For example, recoil events could introduce scatter into the tight correlations between BH 
mass and galaxy bulges [19, 10, 11, 13]. The sudden displacement of the central AGN from these 
merging systems also could affect starburst properties or the gas and stellar distributions [13, 11]. 

We conduct hydrodynamic simulations of galaxy mergers to follow the dynamics and accretion 
of recoiling BHs in evolving merger remnant potentials. Due to the many free parameters involved, 
we undertake a large parameter study with dozens of galaxy merger models and a wide range of 
kick velocities. We are thus able to observe trends in the behavior of recoiling BHs in different 
environments. Our simulation methods and merger models are outlined in § 2. We discuss the 
trajectories of recoiling BHs in § 3. 1. In § 3.2, we examine the total and offset lifetimes of recoiling 
AGN, and in § 3.3, we examine the effects of GW recoil on the co-evolution of SMBHs and their 
host galaxies. We conclude in § 4. 

2. Methods 

We simulate galaxy mergers using the smoothed particle hydrodynamics (SPH) code GADGET- 
3 [20]. The code includes models for radiative cooling, star formation, and supernova feedback 
[21], as well as BH accretion, as described in [22]. In order to model recoiling BHs in these 
simulations, we allow for an arbitrary kick velocity Vk to be aded to the remnant BH at the time 
of the BH merger, fmrg- In most of our analysis, this velocity is scaled to the central escape speed 
at the time of the merger, v^cit^g). We allow the BH to accrete both from ambient gas via the 
Bondi-Hoyle formula and from an ejected disk of gas with a time-dependent accretion rate. These 
prescriptions are described in more detail in [11]. 

We have tested a total of 62 different galaxy merger models in which we vary the galaxy mass 
ratio (q), the total galaxy mass, the gas fraction (/gas), and the orbital configuration. For easier ref- 
erence to these models throughout the text, we assign each a name given by q[value]fg[value][orb], 
where "q" and "fg" denote the galaxy mass ratio and initial gas fraction, and each letter orb is iden- 
tified with a specific orbital configuration. (High- and low-total-mass models are not discussed 
explicitly here; see [11].) For each model, we simulate both a merger with no recoil kick and a 
merger with Vk/v esc = 0.9. For a subset of these models, we also simulate intermediate values of Vk- 
We refer the reader to [1 1] for the full details of our initial conditions. 
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Figure 1: Comparison of recoil trajectories in gas-rich (left panel, model q0.5fg0.6a) and gas-poor (right 
panel, model q0.5fg0.1a) mergers. In each case, the recoiling BH separation from the galactic center is 
plotted versus time for varying v'k/v esc within a single galaxy model. Color-coded values of Vk/v e sc are 
indicated on the plots. The x-axis is the time after the BH merger, t — fmrg- 



3. Results 

3.1 Recoiling BH Trajectories 

Fig. 1 demonstrates that the gas content of galaxies greatly influences recoiling BH trajectories. 
In the left panel, trajectories with Vk = 0.4 — 1 .0 v esc are shown for a merger in which the progenitor 
galaxies had 60% gas initially, while the right panel shows the same set of Vk/v esc trajectories, but 
for galaxies with 10% gas initially. In both cases, BH kicked with 0.8 < Vk < 1.0 v esc are still on 
large orbits by the end of the simulation, 2.9 Gyr after the recoil. However, the BH displacement 
is smaller in the gas-rich case, and BHs with Vk < 0.7 v esc (640 km s _1 in this merger model) are 
confined to the central few kpc of the galaxy. Thus, in gas-rich mergers, even relatively large kicks 
may fail to eject the BH from the central region. 

Recoiling BH dynamics also depend indirectly on gas content via the inflow of cold gas into the 
center of the merger remnant. We find that nearly-equal-mass, gas-rich mergers may experience 
a rapid increase in central v esc during the final coalescence of galaxies. Because the BH is also 
expected to merge during this time, the BH dynamics may depend sensitively on the timing of 
these two events. However, the effect is much smaller for mergers with smaller q and / gas . 

3.2 AGN Lifetimes 

We quantify the observability of recoiling BHs by calculating their lifetimes as either spatially- 
or kinematically-offset AGN. We use three different definitions of an AGN: Lboi > (10%LEdd, 
3%LEdd> & 3.3 x 10 9 L ), where Lt, i and L^dd are the bolometric and Eddington luminosities. 
For kinematically (spatially) offset AGN, we require a BH velocity offset vbh > V min (Rbh > Rmin) 
from the stellar center of mass. Figs. 2a & b show AGN lifetimes (?agn) for Vmin = 500 km s _1 and 
^min = 1 kpc, respectively. 

Figs. 2c & d show the evolution of Lt, i, Rbh, and vbh for two examples of recoil (vk/v esc = 
0.7 & 0.9) in the q0.5fg0.6a merger model. These illustrate two distinct physical scenarios for 
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(a) - Kinematically-offset AGN lifetimes (b) - Spatially-offset AGN lifetimes 





Figure 2: Figs. 2a & b show the kinematic ally- and spatially-offset ?agn for v m ; n = 500 km s~'and R m i n = 1 
kpc, respectively. The 4 galaxy models shown are color-coded according to the plot labels. The dashed 
lines and black circles mark the simulations chosen as examples to illustrate the recoiling AGN luminosity 
and dynamics in Figs. 2c & d. Fig. 2c shows Lb i, ^bh. an d 1 'bh f° r the simulation with model q0.5fg0.6a 
and Vfc/vesc = 0.7. The vertical red line marks the time of the BH merger. The horizontal red lines in the Lbol 
plot mark the three AGN definitions described in the text. The horizontal red lines in the 7?bh and vbh plots 
denote R m i n & v m j n , respectively. Fig. 2d shows the same data as Fig. 2c, but for Vk/v e sc = 0.9. 



kinematic offsets. For recoils near v eS c> the BH carries along a small amount of gas, its luminosity 
decreasing monotonically with time. For intermediate Vk/v eS c> the BH may make repeated passages 
through the central gaseous region. In massive galaxies, v esc may be large enough for BHs on these 
bound orbits to exceed v m i n while renewing their fuel supply. Both mechanisms for producing 
kinematically-offset AGN result in ?agn up to a few 10s of Myr in the merger models shown. 

Spatially-offset AGN with /?bh > 1 kpc (Fig. 2b) correspond to an angular separation of ~ 0.55 
arcsec at z = 0. 1 ; these could be resolved with the HST and JWST if not observed edge-on. In some 
cases, spatially-offset AGN (with /?bh > 1 kpc) may have lifetimes of up to ~ 100 Myr, but they 
are unlikely to appear as bright AGN. Figs. 2c & d illustrate why this is so. In Fig. 2c, the largest 
BH separations correspond to minima in L\, i. In Fig. 2d, the BH displacement is much larger, but 
the BH has no fuel supply for most of its post-merger evolution. 

Finally, we note that in addition to producing offset AGN, GW recoil alters the total AGN 
lifetime. Figs. 2c & d show that the bright AGN phase is cut short by the recoil event, but in the 
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Figure 3: Mbh — 0* relation for simulations with no recoil kicks (black circles) and with Vk/v esc = 0.9 (red 
triangles). Mbh is the BH mass at the end of the simulation (fmrg + 2.9 Gyr) and (7* is the stellar velocity 
dispersion averaged over 100 random sight lines, where error bars give the range of sampled values. The 
solid black and red dashed lines are least-squares fits to the no-recoil and high-recoil data, respectively. The 
fit parameters are indicated on the plot. 

Vk/vesc = 0-7 simulation, accretion continues a lower levels for an extended time. In fact, the low- 
luminosity ?agn in this case is actually longer than for an equivalent stationary BH. This result is 
generic to many of our intermediate-v^ simulations. 

3.3 Co-evolution of BHs and Galaxies 

Because recoiling BHs have shorter lifetimes as bright AGN, they generically have smaller 
final masses than stationary BHs, in some cases by up to a factor of ~ 5. In Fig. 3, we compare 
the Mbh — relation resulting from our sets of simulations with no recoil and with large recoil 
(vk/vesc = ^ e stress tnat tne data s h° wn in Fig. 3 are not expected to reproduce the observed 
Mbh — <7* relation [23, e.g.], as we have varied parameters systematically. In particular, galaxies 
with smaller v esc (and smaller a*) will more frequently experience large Vk/v esc recoils, which 
could steepen the Mbh — slope relative to that shown here. However, by comparing the relative 
differences between the correlations in Fig. 3, we gain insight into how individual GW recoil events 
may contribute to BH/galaxy co-evolution. 

Fig. 3 shows that GW recoil lowers the offset of the Mbh — c* relation and increases scatter. 
While the former effect is somewhat model-dependent, an increase in intrinsic scatter is a necessary 
consequence of recoil. This is because GW recoil introduces another degree of complexity into 
the determination of the final BH mass. The scatter in the Vk/v esc = 0.9 population in Fig. 3 is 
almost a factor of two larger than in the = population (0.24 and 0.13 dex, respectively). In a 
cosmological framework, scatter would increase even further as some recoiling BHs were replaced 
by new BHs via subsequent mergers. 

In addition, GW recoil can directly affect the central structure of galaxy merger remnants. The 
left panel of Fig. 4 compares the star formation rates (SFRs) of simulations with v^/v esc = & 0.9, 
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Figure 4: Left panel: Total SFR vs. time for the qlfg0.4a merger model. The simulations shown have v\ = 
(black dashed line), v^/v esc — 0.9 (red solid line), or no BHs at all (blue dotted line). The vertical line marks 
the time of BH merger in the simulations with BHs. Right panel: Final stellar density profiles for the same 
simulations, with the same color coding. The thick curves are the "new stars" that formed during the merger, 
and the thin curves are the "old stars" from the progenitor galaxy disks. Poisson error bars are shown for each 
radius bin. The vertical dotted line denotes the gravitational softening length (80 pc) in these simulations. 

and one with no BHs. The simulation with large recoil clearly has a higher post-merger SFR than 
the no-recoil simulation, similar to the case with no BHs. This is a direct result of the displacement 
of AGN feedback energy via GW recoil, which leaves a supply of cold gas in the central region. 
Owing to this enhanced SFR, a steeper density cusp forms in the merger remnant. The right panel 
of Fig. 4 shows the final stellar density profiles in these remnants; the newly-formed stars dominate 
within the central kpc. The simulations with recoil and without BHs have about twice the central 
density of the stationary-BH case, which corresponds to a 3% increase in total stellar mass. 

4. Conclusions 

We have conducted the first detailed parameter study of GW recoil in hydrodynamic simula- 
tions of galaxy mergers. Our large suite of simulations allows us to analyze trends in GW recoil 
events with varying galaxy merger parameters and recoil velocities. We estimate the observable 
lifetimes of recoiling AGN, as well as the impact of GW recoil on BH/galaxy co-evolution. 

We find that galaxy gas content and central escape speed at the time of the merger are the most 
important factors in determining the trajectory of a recoiling BH. However, the latter may vary 
rapidly during the time of BH coalescence, especially in gas-rich, nearly-equal-mass mergers. This 
adds an element of unpredictability to the recoiling BH dynamics. 

Recoiling AGN may be distinguished from their stationary counterparts via kinematic or spa- 
tial offsets. In our models, recoiling AGN with vbh > 500 km s _1 have ?agn up to a few tens of 
Myr; those with /?bh > 1 kpc have ?agn up to ~ 100 Myr, though generally at lower luminosities. 
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GW recoil also reduces the total lifetimes of bright, merger-triggered AGN, but may increase the 
low-luminosity ?agn if the BHs resume accretion upon settling back to the center. 

GW recoil may disrupt the coordinated growth of galaxies and BHs. Accordingly, the Mbh — 
a* relation resulting from our simulations with rapidly recoiling BHs has almost twice as much 
intrinsic scatter as that of an equivalent sample with stationary BHs. This additional scatter is an 
unavoidable outcome of GW recoil events. We have also shown that recoiling BHs may affect their 
host galaxies directly, owing to the displacement of AGN feedback from the galactic center. This 
extends the starburst phase and results in a merger remnant with a steeper central cusp. 

In conclusion, our results indicate that GW recoil may be a non-negligible addition to the 
standard picture of BH/galaxy co-evolution. Further, the substantial lifetimes derived for some 
offset AGN suggest that observing such objects might be a real possibility in the near future. 
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